New class of self-similar solutions for vacuum plasma expansion 
admitting monoenergetic ion spectra 
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We report a new class of self-similar solutions for plasma expanding into vacuum that allows 
for quasi-monoenergetic ion spectra. A simple analytical model takes into account externally con- 
trolled time-dependent temperature of the hot electrons. When the laser temporal profile is tailored 
properly, the quasi-neutral self-similar expansion of the plasma results in ion concentration in the 
phase-space at a particular velocity thus producing a quasi-monoenergetic spectrum. We prove this 
analytical prediction using a ID partice-in-cell (PIC) simulation where the time-dependent plasma 
temperature is controlled by two laser pulses shot at a foil with a suitable time delay between them. 
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Generation of highly energetic ions and protons beams 
from the laser interaction with thin foil targets in rela- 
tivistic regime (/ > 10^* W cm^^ fim^, where / is the 
intensity of the laser and A its wavelength) is one of the 
fast developing research fields [l|, [3, [3, 0, H, 0| . Particu- 
larly interesting is the generation of quasi-monoenergetic 
ion beams, because these have a number of important 
potential applications in medical physics, inertial con- 
finement fusion, compact ion accelerators ^7i.,8i.,9i.il0.]. 

When the laser pulse strikes a thin solid density foil, it 
heats electrons at the front surface of the target. These 
hot electrons traverse the target and leave it at the 
rear side. A space charge is build up that leads to a 
huge electrostatic potential. This, in turn, accelerates 
the background ions to high energies. This mechanism 
is known as target normal sheath acceleration (TNSA) 
d, i, H, i, 0, MM E^- From theoretical point of view 
the generation of energetic ions from thin foil targets has 
been studied well in past based on the approach of self- 
similar expansion of plasma 11 , [l3, [3] ■ In this model, 
the plasma fluid is quasi-neutral, ions and electrons both 
have similar flow on the ion sound time scale. In the 



works [14|, [15| also kinetic effects at the front of the ex- 
panding plasma were included in this self-similar model 
for different geometries. 

All the self-similar models up to now have predicted 
broad band ion spectra, decaying exponentially at high 
energies. Yet, most of practical applications require mo- 
noenergetic ion beams. A straightforward way to over- 
come this difficulty is to select a group of ions with the 
same energy, e.g. using ions with the same initial condi- 
tions. Indeed, when a monolayer of protons, or a proton- 
rich dot is placed on the back surface of a higher-Z tar- 
get, then all the light protons are accelerated nearly to 
the same energy [3, [^, because all the protons in the 
monolayer feel same electric field. Still, the number of 
accelerated protons out of the monolayer and thus the ef- 



ficiency of the process are very limited. Recently, it was 
also predicted that ion acceleration from multi-species 
targets may result in peaked spectra 0, [13] ■ 

In this Letter, we report a new type of self-similar solu- 
tion for the quasi- neutral vacuum-plasma expansion that 
does allow for monoenergetic spikes in the ion spectrum. 
As we will see, this new self-similar solution requires a 
suitable time-dependent temperature of laser-generated 
hot electrons. Experimentally, this corresponds to a tai- 
lored temporal profile of the driving laser pulse, because 
the temperature of the hot electrons is a definite function 
of the laser intensity 1^, 17 1. The rising electron tem- 
perature with time is accompanied by a rising electro- 
static field inside the expanding plasma. Consequently, 
this results in a higher acceleration rate of ions at later 
times. These rear ions catch up with the leading ions, 
which began to accelerate earlier. When this catching 
up occurs, a bunch of ions concentrated in phase space is 
formed. If driven too harsch, this catching up may result 
in "wave breaking/' [l8| . However, when the temperature 
is changed gently, one may stay "just at the verge" of 
wavebreaking. This is exactly the solution we present 
here. 

Under the quasi- neutral assumption i.e Ug « Zui 
where Ug and Ui are the electron and ion densities repec- 
tively, and Z is the ionization state, the dynamics of the 
plasma expansion is governed by the continuity equation 
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and the Euler equation for the ion fluid 
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Here Vi is the ion fluid velocity, is the ion mass, 
T^(t) = To(t)(ne(0,t)/ne(0,0))'''-i is the local electron 
temperature and 7 is the expansion polytrope index for 
electrons. We also assume that the plasma electrons obey 
the Boltzmann statistics 
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ne{r,t) = ne{t)exp{e4>/Te 



(3) 



where ne(t) is the time-dependent electron density, — e is 
the electron charge, and 4>{r, t) is the electrostatic poten- 
tial. The ions are assumed to be cold initially. 

It is well known that the above system of equations 
permits the following self-similar solution 



Me) 



R^, ^ = x/R, 



n, « Zn, = nooiR.o/R) N,{0, N,{0) = 1, 



(4) 



where R{t) is the time dependent charactersitic size of 
the plasma, ^ is the similarity coordinate {x is the posi- 
tion coordinate), and Ni is a positive unknown function. 
The subscripts i and e refer to the ions and electrons 
respectively. The rationale behind this self-similar solu- 
tion is that a linear velocity-radius relation is the correct 
limit for the asymptotic stage of the expansion, whenever 
the R(t) greatly exceeds the initial value Rq = R{0) fl^ . 
We may see that this solution rules out the possibility of 
generation of quasi-monoenergetic ion beams because the 
ion velocity scales linearly with the similarity coordinate 

e- 

In this Letter, we introduce a more general self- 
similarity ansatz and assume that the plasma expansion 
scales as 



that self-similar solutions are not well defined at the ori- 
gin and must be considered at some distance from the 
origin only 

One easily retrieves the old similarity ansatz ([4]) by 
setting C = in ([6]). Thus, all the previous self-similar 
solutions for plasma expansion into vacuum were missing 
the second term in Yet, this term is very important. 
It is the second term in ^ with the non-zero constant 
C that allows for a new class of self-similar solutions and 
gives the option for quasi-monoenergetic ion spectra. 

The condition to have a monoenergetic spike in the 
ion spectrum is the presence of a stationary point in the 
expansion velocity, dvi{(_-ai) / dS^ ~ 0, that gives 



C 



d^ 



(7) 



The condition ([7]) means that there is a group of ions 
around the point which move all with the same veloc- 
ity 



Rf (Cn 



(8) 



This ion beam produces a spike in the energy spectrum 
around — miv'^/2. 

Now we substitute the formula ^ in Eq. ^ and obtain 



m ^ Nf[N,(-C{l + a)] 
d^ [aC^ - ' 

where we have used the following relations 



(9) 



Tie ~ Zn, = noo{Ro/R) N,{0, N,{0) = 1, (5) 

where /(^) is a function which is yet to be determined. 
Because the new solution must satisfy the equation of 
continuity, we have from Eq.([T]) 



brriiRo V R 



(10) 



a and 6 are constants. It follows from Eqs. ^ and ^ 
that the density of the monoenergetic ion beam is (in the 
isothermal case, 7 = 1) 



/(0-e-c/iv,(e). 



(6) 



Here C is a constant of integration. Physically, the fluid 
velocity Vi must be monotonic. Consequently, the so- 
lution ^ has a meaning only in the range ^ < ^max 
where ^max marks the front end of the plasma expan- 
sion in vacuum. We require that C < ^max^^i (Cmax)- 
This restriction means that C is small, because Ni decays 
exponentially with ^. Thus, the new solution is hardly 
distinguishable from the classic one at small ^. Only ap- 
proaching ^max the solution is modified significantly. For 
the sake of accuracy, we mention that our solution leads 
to a non-zero fluid velocity at the origin, Vi{0) = —CR. 
This velocity, however, is exponentially small (just as C) 
and has no great physical meaning. It is well known, 



(em) = {C/2b) U + ~ 46) 



,1/2' 



and its velocity 



R(u-2b (em + (em"4fe) 



-1/2 



(11) 



(12) 



It is clear that the density and energy of the monoen- 
ergetic bunch do not depend on the constant a. Hence, 
in principle this constant can be chosen freely. The con- 
stant b can be expressed in terms of the constant C and 
density of the monoenergetic bunch, which is sought ex- 
perimentally, as 



(13) 
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One can draw further conclusions regarding ion accel- 
eration from the self-similar solution. We see from the 
Eq. that the ion acceleration is inversely proportional 
to i?, which is an increasing function of time. We obtain 
the expression for R on integrating the first relation in 
Eq.dini) 



R 

Rq 



a Rq , 



(14) 



where Vj: = u/ (i?o/-R/)^^"- Here Vf and Rf are the 
velocity and size of the expanding plasma, in the case of 
a time independent electron temperature, respectively. 
These can be determined from the results of Ref.[l3|. 
Once the expression for R is obtained, we can determine 
the required electron temperature profile from the second 
part of the Eq. pO)) : 
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where Tf = rriiV^ is the time independent electron tem- 
perature. This relation describes the electron tempera- 
ture variation in terms of the constants a and b. These 
constants define how time dependent electron tempera- 
ture profile must be tailored to produce the monoener- 
getic ion beam. The idea of controlling the plasma elec- 
tron temperature has been recently studied by means of 
Vlasov and PIC simulations [3] ■ It is found that a suit- 
able temperature variation of the plasma electrons indeed 
can produce the mono-energetic ions. 

An analytical solution of Eq. Q is hardly possible and 
we integrate it numerically. Fig[T] shows Ni and Vi as a 
function of ^ for the parameter set a = 5, 6 = 0.19, C = 
3.3636 X 10^*^,and7 — 1. FiglT] describes the self-similar 
expansion of the quasi-neutral plasma. The most im- 
portant result is the behaviour of Vi showing flattening 
with ^ near ^ = 2. This flattening manifests the produc- 
tion of monoenergetic ions. The velocity of expansion 
reaches its maximum at the flattening point. It corre- 
sponds also to the position of the ion front, where the 
self-similar solution has to be terminated. The energy 
spectrum of ions produced by this self-similar solution 
is given in Figl^l We observe the quasi-monoenergetic 
bunch formation from the figure. 

We have performed a ID particle-in-cell (PIC) simula- 
tion of the quasi-monoenergetic protons generation from 
a thin foil target with the help of virtual laser plasma 
laboratory (VLPL) code |20j|. In the simulation, we take 
a thin foil of thickness 1.5A with an initial density ramp 
and solid density 30nc, where Uc is the critical density 
and A is the wavelength of the laser pulse. 

To tailor the time dependent temperature of the hot 
electrons, we use two driving laser pulses shot at the tar- 
get with a time delay between them. The simulation 



FIG. 1: Variation of Ni andwi with ^ showing the self-similar 
flow and velocity of the expanding plasma. The inset is a 
zoomed in view of the velocity of expansion near the flattening 
point. The parameters are given in the text. 
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FIG. 2: Ion energy spectrum corresponding to the self-similar 
solution given by Eq.® for 7 = 1. Here the energy is nor- 
malized by miB? /2 and the ion density Ni is normalized by 
the value at the center of the thin foil target. The other pa- 
rameters are same as in the Fig[T] 



parameters for the laser pulses are; ai — eAi/mc? = 
1, cTi — 5A, a2 — eAijnK? = 3, cti = 15A and A = 
0.8 /im, where t\ and T2 are the pulse durations, and 
a\.2 are the normalized vector potentials of the two laser 
pulses. Here subscripts 1 and 2 represent the first and 
second laser pulse respectively. The time delay between 
the two pulses is 60 laser periods. Both laser are y- 
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FIG. 3: The phase space evolution of protons at the various 
stages of thin foil expansion. The middle and lower panels of 
the figure describe the momenta flattening, which ultimately 
leads to the generation of monoenergetic protons from the 
thin foil target. 



polarized and move in x direction. The resolution in x- 
direction is O.OIA while 200 numerical particles per cells 
are chosen to run the simulation. 

Fig [3] depicts the phase space evolution of protons at 
the various stages of plasma expansion. The upper panel 
depicts the stage of the thin foil expansion under the 
influence of the first laser pulse. It is clear that the sin- 
gle laser pulse leads to the standard self-similar plasma 
expansion given by Eq.(l4|) and thus clearly rules out the 
possibility of concentration of protons in the phase-space. 
The proton velocity scales linearly with the distance from 
the target. Later (at i = 80 • 27r/ti; where w is the carrier 
frequency of the laser pulse), when the second pulse hits 
the target (c.f. middle panel of the Figl3]), the dynamics 
of the expanding plasma changes and one can see mo- 
menta flattening in the phase space. The momenta flat- 
tening forms the monoenergetic protons bunch. Phys- 
ically, the second laser pulse increases the hot electron 
temperature that leads to a rising electrostatic field and 
a larger acceleration of ions at later times. These ions, 
which form a second group, catch up with the first group 
of ions, which began to accelerate earlier. This process 
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FIG. 4: The energy spectrum of the protons showing the 
monoenergetic peak around 4 MeV from the PIC simulation. 
The proton density Ni is normalized by the critical density. 



flattens the velocity profile and forms a bunch of quasi- 
monoenergetic ions . As the time passes, this bunch of 
ions accelerates and acquires substantial energy as shown 
in the lower panel of the FiglJ] Fig[3]shows the ion energy 
spectrum obtained in the PIC simulation. It can be seen 
that there is a significant number of ions with the peak 
energy around 4 MeV and the energy spread is a few per 
cent. The number density of quasi-monoenergetic ions 
can be further enhanced by using higher density targets 
and longer laser pulses. Using the two-pulse technique - 
or, more general, a tailored time-dependent laser inten- 
sity - it will be possible to produce light as well as heavy 
ions {Z > 1) with monoenergetic spectra. However, the 
heavy ion case may demand higher laser intensities. 

In summary, we have found a new class of self-similar 
solutions for the quasi-neutral vacuum-plasma expansion 
that allows for quasi-monoenergetic ion spectra. The 
quasi-monoenergetic spectrum corresponds to the flat- 
tening of the ion momentum dependence on the expan- 
sion coordinate. Thus, ions pile up at a particular veloc- 
ity in the phase space and form a quasi-monoenergetic 
bunch. This new scheme demands for laser pulses with 
a tailored temporal profile. To prove realizability of the 
new scheme, we have performed a ID PIC simulation us- 
ing the code VLPL. In the PIC simulation, we use two 
laser pulses hitting the foil target with a sutiable time 
delay between them. Although the developed model is 
simple, it catches the essence of the physical process quite 
well. The proposed scheme works already for moderate 
laser powers which should make the laser-plasma gener- 
ation of monoenergetic protons (ions) a stable and well 
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controlled process. 
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